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■ Roy's equations are used to check if the scalar-isoscalar tttt scattering ampli- 

tudes fitted to experimental data fulfill crossing symmetry conditions. It is shown 
that the amplitudes describing the "down-flat" phase shift solution satisfy cross- 
ing symmetry below 1 GeV while the amplitudes fitted to the "up-flat" data do 
not. In this way the long standing "up-down" ambiguity in the phenomenological 
determination of the scalar-isoscalar tttt amplitudes has been resolved confirming 
the independent result of the recent joint analysis of the tt + tt~ and data. 
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Abstract 



^ 1 Introduction 

H , 

In 1997 a new analysis of the tc — > n + n n reaction on a polarized target was 
performed in the m 7Tn effective mass range from 600 to 1600 MeV [1]. For the first 
time the pseudoscalar (7r-exchange) amplitude was separated from the pseudovector 
(ai-exchange) amplitude in the region of the the four-momentum transfer squared 
from —0.005 to —0.2 (GeV/c) 2 . Below 1000 MeV, where the S- and P-waves strongly 
interfere, the partial wave analysis of the 7r + 7r~ data provided us with two scalar- 
isoscalar solutions, called "up" and "down", which differ by their intensities. Lack 
of information on a sign difference between the phases of the S- and P-waves near 
the position of the p resonance led us to other two branches of the "up" and "down" 
amplitudes named "steep" and "flat". It was shown in [2] that both "up-steep" and 
"down-steep" S'-wave isoscalar amplitudes significantly violate unitarity below 1 GeV 
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and should be rejected as nonphysical. Two remaining "flat" amplitudes survived 
the unitarity check and other tests were needed to resolve the existing "up-down" 
ambiguity. 

In 2001 new experimental data on the 7r°7r° production from the E852 collaboration 
appeared [3] and were used in a joint analysis of the 7r + 7r~ and 7r°7r° data [4]. The 
7r°7r° data are very useful to compare with the 7r + 7r~ data due to an absence of the 
P-wave in the 7r°7r° channel and therefore a lack of the "up-down" ambiguity. The 
one-pion and ai-exchange model described in [4] was used to calculate the S'-wave 
intensities of the 7r°7r° production by choosing as an input the "up-flat" or "down-flat" 
phase shifts. The isospin relations between the 7t + tt~ — > 7t + tt~ and 7r + 7r~ — > 7r°7r° 
amplitudes supplemented by the parameterization of the isotensor scalar amplitude 
taken from [1] were helpful in these calculations. It was shown that the 7r°7r° S'-wave 
intensities determined for the "down-flat" phase shifts agree with the experimental 
values within the errors. However, for the "up-flat" phase shifts in the m nn range from 
850 to 970 MeV important differences between the calculated 7r°7r° intensities and the 
corresponding experimental values occur. This fact led the authors to a conclusion 
that the "up-flat" data set should also be rejected. 



2 Roy's equations as a test for the nn amplitudes 

Another independent test of the "up-flat" and "down-flat" amplitudes consists in 
checking if they fulfill crossing symmetry conditions below 1 GeV. In order to achieve 
this task we have used Roy's equations [5] for the scalar-isoscalar, £ = 0, / = 0, scalar- 
isotensor, £ = 0, 1 = 2, and the vector-isovector, £ = 1, I — 1, tttt partial waves 
determined in a wide m nw range. We were especially interested in the m nn region 
between 800 and 1000 MeV where differences between phase shifts of the "up-flat" and 
"down-flat" data sets are largest and reach about 45° (see Fig. 4 in [2]). In a recent 
analysis of Roy's equations [6] a special attention was put on the effective mass lower 
than 800 MeV. 

As an input in Roy's equations we have used imaginary parts of the partial waves 
amplitudes //(s) related to the nn phase shifts and inelasticities r\\: 



where /x is the charged pion mass and s = m^ n . 

Below 970 MeV the following Pade representation of both the "up-flat" and "down- 
flat" phase shifts has been taken: 

tan x0 ( „) _ Eto a2i+ik 2l+1 , , 

^W)- U 3 =i{k 2 /a2 ._ iy W 

where k = |v ' s ~~ 4/i 2 is the pion momentum and otj (j = 1, . . . , 7, 9) are constant 
parameters. Above 970 MeV up to 2 GeV our coupled channel model [7] amplitude A, 
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fitted to the "down-flat" data, and the amplitude C, constrained by the "up-flat" data, 
were applied. In the fits we also used the near threshold phase shifts calculated from 
the differences Sq — 5\ obtained in the high statistics K ei decay experiment [8]. The 
scattering length a[j and the slope parameter 6° are directly related to the constants 
otj\ a[j = — and = —a\[i (0.5/i~ 2 + a^ 1 + 1 + ctg 1 — afj — a^/z. 

The parameterization of isotensor wave using a rank-two separable potential model 
has been described in [9] where detailed analysis of the present study is presented. 

For the P-wave, from the irir threshold till 970 MeV, we have used an extended 
Schenk parameterization [6]: 

t an Sl(s) = ]jl-^f k 2 (A + Bk 2 + Ck A + Dk 6 ) ( 4 f_~/ P ) , (3) 

where A is the P-wave scattering length and s p is equal to the p-mass squared. Above 
970 MeV we took the i^-matrix parameterization of Hyams et al. [10]. The parameters 
C and D were chosen to join smoothly the phase shifts given by both parameterizations 
around 970 MeV. 

The contributions to Roy's equations from high energies (m n7T > 2 GeV) and from 
higher partial waves (/ > 1) are called driving terms. They are composed of contribu- 
tions from the / 2 (1270) and p 3 (1690) resonances and from the Regge amplitudes for 
the Pomeron, p- and /-exchanges. The Breit-Wigner parameterization with masses, 
widths and tttt branching ratios taken from [11] were used for / 2 (1270) and p 3 (1690). 
For the Regge parts we have used formulae of [6] without the ^-crossed terms. We have 
found that the / 2 (1270) resonance dominates in the scalar isoscalar wave and that the 
introduction of the p3(1690) has a significant influence on the isotensor and isovector 
waves. In the isoscalar wave the Regge contributions are more than 10 times smaller 
than the resonance contributions but for the isospin 1 and 2 they are of the same order. 

The thirteen constants, six for the scalar-isoscalar wave in (2), four for the isotensor 
wave and three for the isovector wave in (3), were calculated from the simultaneous 
fits to data and to Roy's equations separately for the "up-flat" and "down-flat" data. 
We have used the CERN MINUIT program with the \ 2 test function defined by 
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1 2 



, (4) 



where ipj (s,) and Atpj (sj) represent the experimental phase shifts and their errors, 
respectively, Sj = [4j + 0.001]/i 2 for j = 1, 11 and si 2 = 46.001/i 2 . The real parts 
Re f\ n have been calculated from (1) under an assumption that the inelasticity r\\ is 
equal to 1 and the phase shifts 5f are equal to (pj(sj). Other real parts, denoted by 
Re fo Ut , constitute the output values calculated from Roy's equations. We take a Af 
value of 0.5 x 10~ 2 to obtain acceptable fits to Roy's equations. 18 experimental values 
of the "up-flat" or "down-flat" data between 600 and 950 MeV were used in addition 
to six data taken from [8]. 
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In Fig. la and lb we present results of fits to the "up-flat" and "down-flat" phase 
shifts and to Roy's equations (solid lines). In both cases differences | Re fl ut - Re f( n \ 
were of the order of 10~ 3 in all three partial waves. The x 2 fc> r 18 points between 600 
and 970 MeV was 16.6 in the "down-flat" case and as large as 46.4 in the "up-flat" 
one. We see in Fig. la that the solid line lies distinctly below the "up-flat" data points 
between 800 and 970 MeV. In contrary, the corresponding line for the "down-flat" case 
in Fig. lb is very close to experimental data in the same range of m nn . In order to 
improve a fit to the "up-flat" data we have used constraints given by the good fit to the 
"down-flat" data. Two parameters were fixed by choosing the values of the scattering 
length and the slope parameter and two others by the values of phase shifts calculated 
from this fit at 500 and 550 MeV. A new fit with these constrains gave an improved 
value of x 2 — 13 for 18 "up-flat" data points, corresponding to the first part of x 2 
in (4) but provided us with an enormous value of x 2 — 1-2 x 10 4 for the second part 
related to Roy's equations. The phase shifts for this amplitude are presented in Fig. 
la by the dotted line. It is clear that a simultaneous good fit to the "up-flat" data and 
to Roy's equations is impossible. 

Apart of the fits to the "up-flat" and "down-flat" experimental points we have also per- 
formed fits to points shifted upwards and downwards by their errors. In these fits the 
same four constraints described above were used below 600 MeV. Up to 937 MeV in 
the "down-flat" case in Fig. Id the curves labeled higher "in" and lower "in" form a 
band including inside a band delimited by the lines higher "out" and lower "out". All 
the curves lying inside these bands correspond to the amplitudes fulfilling the crossing 
symmetry so the "down-flat" data can be accepted as physical ones. In the "up-flat" 
case in Fig. lc the output band lies outside of the input band from 840 to 970 MeV. 
It means that in this case the crossing symmetry is violated by the amplitudes fitted 
to the "up-flat" data. 

In Fig. 2 we have presented the output results for the isotensor and isovector waves 
in the "down-flat" case only since in the "up-flat" case the curves are very similar. The 
"in" curves were not plotted because they are almost indistinguishable from the "out" 
ones. 

3 Conclusions 

We have used Roy's equations as a tool to test if the amplitudes fitted to the 
"up-flat" and "down-flat" phase shifts extracted from the n~p^ — > -n +r n~n data fulfill 
crossing symmetry conditions. We have found that only the S'-wave isoscalar amplitude 
corresponding to the "down-flat" data set can be accepted. The amplitude constrained 
to the "up-flat" data does not satisfy Roy's equations and should be rejected as non- 
physical. This conclusion is in agreement with the independent results obtained from a 
joint analysis of the n + 7i~ and the 7r°7r° production data [4]. In this way a long stand- 
ing "up-down" ambiguity in the nn experimental data has been resolved in favour of 
the "down-flat" data set. 
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Figure 1: a) and c) correspond to the "up-flat" case, b) and d) correspond to the 
"down-flat" case. Fits to the scalar-isoscalar phase shifts of [1] and to Roy's equations 
are denoted by solid lines in a) and b) . Dotted line in a) between two dashed lines rep- 
resents fit with constraints described in the text. Dashed lines in a) and b) represent 
fits to phase shifts moved upwards and downwards by their errors; the corresponding 
lines in c) and d) are called higher and lower, respectively. Lines in c) and d) cor- 
respond to real parts of input amplitudes ("in") and real parts calculated from Roy's 
equations ("out"), all multiplied by 2ks~ 1 ^ 2 . Diamonds represent the K e4 data [8]. 
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Figure 2: Real parts of isotensor a) and isovector b) tttt amplitudes (multiplied by 
2ks~ 1 / 2 ) fitted to the "down- flat" data. Triangles in a) denote data of [12]. Crosses in 
b) are the pseudo-data calculated from the if-matrix fit of [10]. 
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